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We study the efTect of population inversion associated with the electron and hole injection in 
graphene p-i-n structures at the room and slightly lower temperatures. It is assumed that the 
recombination and energy relaxation of electrons and holes is associated primarily with the interband 
and intraband processes assisted by optical phonons. The dependences of the electron-hole and 
optical phonon effective temperatures on the applied voltage, the current-voltage characteristics, 
and the frequency-dependent dynamic conductivity are calculated. In particular, we demonstrate 
that at low and moderate voltages the injection can lead to a pronounced cooling of the electron-hole 
plasma in the device i-section to the temperatures below the lattice temperature. However at higher 
voltages, the voltage dependences can be ambiguous exhibiting the S-shape. It is shown that the 
frequency-dependent dynamic conductivity can be negative in the terahertz range of freqiencies at 
certain values of the applied voltage. The electron-hole plasma cooling substantially reinforces the 
effect of negative dynamic conductivity and promotes the realization of terahertz lasing. On the 
other hand, the heating of optical phonon system can also be crucial affecting the realization of 
negative dynamic conductivity and terahertz lasing at the room temperatures. 



I. INTRODUCTION 

The gapless energy spectrum of electrons and holes in 
graphene layers (GLs), graphene bilayers (GBLs), and 
non-Bernal stacked multiple graphene layers (MGLs) 

opens up prospects of creating terahertz (THz) lasers 
based on these graphene structures. In such structures, 
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GLs and MGLs with optical (3-[l0| and injection 
pumping can exhibit the interband population inversion 
and negative dynamic conductivity in the THz range of 
frequencies and, hence, can serve as active media in THz 
lasers. The most direct way to create the interband pop- 
ulation in GLs and MGLs is to use optical pumping [4] 
with the photon energy hflo corresponding to middle- 
and near- infrared (IR) ranges. In this case, the electrons 
and holes, photogenerated with the kinetic energy hno/2, 
transfer their energy to optical phonons and concentrate 
in the states near the Dirac point [1, [H, [l^ . The amplifi- 
cation of THz radiation fror n op tically-pumped GL struc- 
tures observed recently [13 is attributed to the inter- 
band stimulated emission. However, the optical pumping 
with relatively high photon energies exhibits drawbacks. 
First of all, the optical pumping, which requires complex 
setups, might be inconvenient method in different appli- 
cations of the prospective graphene THz lasers. Second, 
the excessive energy being received by the photogener- 
ated electro-hole plasma from pumping source can lead 
to its marked heating because of the redistribution of the 
initial electron and hole energy MIq/2 among all carries 
due to rather effective inter-carrier collisions. The latter 
results in a decrease of the ratio of the quasi-Fermi ener- 
gies /ie and /i/i to the electron-hole effective temperature 
T that, in turn, complicates achieving of sufficiently large 
values of the dynamic conductivity. As demonstrated re- 
cently the negative conductivity at the THz fre- 
quencies is very sensitive to the ratio of the photon en- 
ergy Ml and the optical phonon energy fiojo , as well as to 



the relative efficiency of the inter-carrier scattering and 
the carrier scattering on optical phonons. The decay of 
nonequilibrium optical phonons also plays an important 
role. 

The abovementioned complications can be eliminated 
in the case of pumping resulting in the generation in GLs 
electrons and holes with relatively low initial energies. 
This in part can be realized in the case of optical pump- 
ing with MIq/2 < fiwo'fl^ ■ Taking into account that in 
GLs huJo ~ 0.2 eV, in the case of CO2 laser as a pumping 
source, fio/^o — 0.5. As shown [l^, in such a case, the 
electron- hole plasma can even be cooled, so that T < Tq, 
where To is the lattice (thermostat) temperature. An- 
other weakly heating or even cooling pumping method 
which can provide low effective temperature T (includ- 
ing T < To) is the injection pumping of electrons from 
n-section and holes from p-section in GL and MGL struc- 
tures with p-i-n junctions. 

In this paper, we study the injection phenomena in 
GL and MGL p-i-n structures and calculate their char- 
acteristics important for THz lasers. The idea to use p-n 
junctions in GLs was put forward and briefly discussed 
by us previously [Tl| . Here we consider more optimal de- 
signs of the structures (with a sufficiently long i-section) 
and account for realistic mechanisms of recombination 
at elevated temperatures (at the room temperature and 
slightly below). 

The paper is organized as follows. In Sec. II, we de- 
scribe the device structures under consideration and prin- 
ciples of their operation. The pertinent equations of the 
model governing the balance of electrons, holes, and op- 
tical phonons (rate equations) are presented in Sec. HI. 
These equations are reduced to an equation governing the 
electron-hole effective temperature. The solution of this 
equation in Sees. IV and V (both analytically in limiting 
cases and numerically) allows us to find the effective tem- 
perature of optical phonons and the current as functions 
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FIG. 1. (Color online) Schematic view of the cross-sections 
of MGL laser structures (a) with chemically doped n- and p- 
sections and (b) with such section electrically induced by the 
side gate-voltages Vp = —Vg and Vn — Vg > 0. 




FIG. 2. (Color online) Band profiles of a GL in p-i-n junctions 
(a) atV = and (b) at forward bias V > 0. Opaque and open 
circles correspond to electrons and holes, respectively. Wavy 
arrows show some interband (recombination) and intraband 
transitions assisted by optical phonons. Smooth arrows indi- 
cate tunneling and thermionic leakage processes. 



of the applied voltage, the structural parameters, and 
the lattice temperature. In Sec. VI, the obtained char- 
acteristics of the injected electron-hole plasma are used 
to calculate the dependence of the dynamic conductivity 
of the latter as a function of the signal THz frequency 
and other quantities. Sec. VII deals with the model lim- 
itations and discussion. In Sec. VIII, we draw the main 
conclusions. 



II. DEVICE MODEL 

We consider devices which comprise a GL or an MGL 
structure with several non-Bernal stacked GLs. It is as- 
sumed that the sections of GLs adjacent to the side con- 
tacts are doped (p- and n-sections) . The device structure 
under consideration is shown in Fig. 1(a). The dc volt- 
age V is applied between the side contacts to provide the 
forward bias of the p-i-n junction. Due to doping of the 
side sections with the acceptor and donor sheet concen- 
tration Ei, the electron and hole Fermi energies counted 
from the Dirac point are fi^ — fJ"i and fih = —fJ-i, where 
fii ~ ft Wi? v^ttS^, h is the reduced Planck constant, and 
Vp ~ 10^ cm/s is the characteristic velocity of the car- 
rier spectrum in GLs. Instead of doping of the side sec- 
tions, the p- and n-sections can be created using highly 
conducting gates over these sections to which the gate 
voltages Vp = —Vg < and Vn — Vg > are applied \U\ . 
In this case, the chemically doped p- and n-sections are 
replaced by the electrically induced sections (see, for in- 
stance, as shown in Fig. 2(b). In single-GL struc- 
tures, Hi cx \/Vg/Wg, where Wg is the spacing between 
the GL and the gate. In the case of MGL structures, the 
situation becomes more complex due to the screening of 
the transverse electric field in GLs although the ef- 
fective electric doping can be achieved in MGL structures 
with about dozen GLs. For definiteness, in the following 
we shall consider primarily the devices with chemically 
doped p-i-n junctions. As shown below, under certain 
conditions the frequency-dependent dynamic conductiv- 
ity of the GL or MGL structures can be negative in a 
certain range of the signal frequencies. In this case, the 
self-excitation of THz modes propagating in the substrate 
serving as a dielectric waveguide (in the direction perpen- 
dicular to the injection current) and lasing is possible. 
The metal gates in the devices with electrically induced 
p-i and i-n junctions can also serve as the slot-line waveg- 
uides for THz waves. Apart from this, the self-excitation 
of surface plasmons (plasmon-polaritons) is possible as 
well (see, for instance, Refs. ^.M Isl-fToj) 

When the p-i-n junction under consideration is forward 
biased by the applied voltage V, the electrons and holes 
are injected to the i-section from the pertinent doped 
side sections. The injected electrons and holes reaching 
the opposite doped section can recombine at it due to 
the interband tunneling or escape the i-section due to 
the thermionic processes. The band profiles in the struc- 
tures under consideration at = and at the forward 
bias V > are shown in Fig. 2. Since the probability 
of such tunneling is a very sharp function of angle of in- 
cidence, the leakage fiux due to the tunneling electrons 
(holes) is much smaller that the fiux of injected electrons 
(holes). The currents associated with the tunneling and 
thermionic leakage of electrons at the p-i-junction and 
holes at the i-n-j unction depend on the electric field at 
the pertinent barriers and the applied voltage V. Because 
the relative role of the leakage currents diminishes with 
increasing length of the i-section 2L, width of the p-i- and 
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i-n- junctions I, and the barrier height at these junctions 
fii, we shell neglect it. The pertinent conditions will be 
discussed in the following. Thus it is assumed that the 
main fractions of the injected electrons and holes recom- 
bine inside the i-section. The recombination of electrons 
and holes in GLs at not too low temperatures is mainly 
determined by the emission of optical photons [l^l . Con- 
sidering the sub-threshold characteristics (i.e., the states 
below the threshold of lasing) and focusing on the rela- 
tively high-temperature operation, we shall account for 
this recombination mechanism and disregard others 
[23 | including the mechanism [25'| associated with the tun- 
neling between the electron- hole puddles (if any) [26l - l29| . 
Due to high net electron and hole densities in MGL struc- 
tures with sufficient number of GLs, the latter mechanism 
can be effectively suppressed [25|. We also assume that 
the net recombination rate in the whole i-section is much 
smaller than the fluxes of injected electrons and holes. 



III. EQUATIONS OF THE MODEL 

Due to rather effective inter-carrier scattering, the elec- 
tron and hole distribution functions (at least at not too 
high energies) can be very close to the Fermi distribution 
functions with quasi-Fermi energies /ig and /i^ and the 
electron-hole effective temperature T. The latter quan- 
tities are generally different from those in equilibrium 
(without pumping) at which ii^ — fJ'h = and T = Tq. At 
the pumping of an intrinsic GL structure, /Ze = — M/i = /i, 
where generally /i > 0. Under these conditions, the quasi- 
Fermi energy /i in the i-section (neglecting the leakage 
and recombination in the lowest approximation) , is given 
by [see Fig. 2(b)] 
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(4) 



""^ are the pertinent characteristic 
times (relatively slow dependent on fi and T), Sq is the 
equilibrium electron and hole density, and Aq is the num- 
ber of optical phonons. Here and in all equations in the 
following, T and Tq are in the energy units. When the 
optical phonon system is close to equilibrium, one can 
put No — [exp{fujJo/To) — 1]""^ = A/'o''. For numerical es- 



timates we set Rq 



Equations (2) and (3) yield the following general formula 
for the structure current-voltage characteristic; 
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2eLI]o 

winter 
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M = eV/2, 



(1) 



Naturally, at 0, T = Tq, so that J\fo = 

[exp{hujQ/TQ) - and J = 0. At V > 0, due to con- 
tributions of the recombination and injection to the en- 
ergy balance of the electro-hole plasma in the i-section, 
the electron-hole effective temperature T can deviate 
from the lattice temperature Tq. The number of opti- 
cal phonons Ao can also be different from its equilibrium 
value A/'o'^. Since Huq is large, in a wide range of temper- 
atures (including the room temperatures) Hcoq ^ Tq 

The electron-hole plasma gives up the energy Hujq 
in each act of the optical phonon emission (interband 
and intraband) and receives the same energy absorb- 
ing an optical phonon. Hence, the net rate of the 
energy transfer from and to the electron hole-plasma 
due to the inreaction with optical phonons is equal to 
2Lfta;o(i?o"*'"'+^o"*™)- Considering Eqs. (3) and (4) and 
taking into account that the Joule power associated with 
the injection current is equal to Q = JV = 2eLR^J^*'^^V, 
an equation governing the energy balance in the electron- 
hole plasma in the i-section can be presented as 



where e is the electron charge. 

The terminal current between the side contacts (per 
unit length in the lateral direction perpendicular to the 
current), which coincides with the recombination current, 
is given by 



J = 2eLR, 







(2) 



The rate of the optical phonon-assisted interband tran- 
sitions (recombination rate) i?™*'^'' and the rate of 
the intraband energy relaxation associated with optical 
phonons i?™**"" can be calculated usin g th e following sim- 
plified formulas [1, (see, also Ref. |20| : 
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Here the left-hand side corresponds to the power received 
by the electron-hole plasma in the i-section from the 
pumping source, whereas the right-hand side correspond 
to the power transferred to or received from the optical 
phonon system. 

The number of optical phonons is governed by an equa- 
tion which describes the balance between their generation 
in the interband and intraband transitions and decay due 
to the anharmonic contributions to the interatomic po- 
tential, leading to the phonon-phonon scattering and in 
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the decay of optical phonons into acoustic phonons. This 
equation can be presented in the form 



(A/'o-A/'o^) 

decay 
"^0 



1 



(TVo + l)exp 



-Mo 



(AAo + l)cxp(-^)-AAo 



(7) 



decay 

can be markedly longer than r™ 



where Tq is the optical phonon decay time. This time 

winter Q^j^d ^^^^^^ 



'0 



, partic- 



ularly in suspended GLs, so that parameter 77Q 



decay 



decay 



j-j-mter ^^^^ excccd Or Substantially exceed unity. 



As shown I30I-I3 



in GLs is in the range of 1 



10 ps. As calculated recently [17|, the charactertic times 
^inter ^mtra longer than 1 ps. If so, the sit- 

uation when tjq^'^"'^ < 1 appears also to be feasible. The 
optical phonon decay time might be fairly short depend- 
ing on the type of the substrate. 

Instead of Eq. (7) one can use the following equation 
which explicitly reflexes the fact that the energy received 
by the electron-hole plasma from the external voltage 
source goes eventually to the optical phonon system: 



decay 

Vo 



eV 



(AAo + l)exp(^ll^)-AAo 



eq 



(8) 

Using Eqs. (6) and (7) or Eqs. (6) and (8), one can find T 
and Ao as functions of V and then calculate the current- 
voltage characteristic invoking Eq. (5), as well as the dy- 
namic characteristics. 
Equation (8) yields 



The ratio 770 



j-intra fainter 
^0 /^O 



;q is actually a function of 
fj, and T. The 770-/1 and r]Q — T dependences are as- 
sociated with the linearity of the density of states in 
GLs as a function of energy. To a good approxima- 
tion these dependences can be described ^16,] by function 
rio = h^ujll{iop? -I- 7r2r2) with r/o cx {huo/Tf at < T 
and 770 oc (fta;o//i)^ at hujf) > /i 3> T. Thus, considering 
Eq. (1), 770 in Eq. (10) is given by 



7?0 



(11) 



Introducing the effective temperature of the optical 
phonon system <d such that Ao = [exp(fta;o/6) ~ l]~^j 



e 



ln(l+AAo-i)' 



(12) 



and substituting Ao from Eq. (9) to Eq. (11), one can 
relate Q and T. Then calculating the T—V dependences 
using Eq. (10), one can find the pertinent Q — V depen- 
dences. 



IV. EFFECTIVE TEMPERATURES AND 
CURRENT- VOLTAGE CHARACTERISTICS 
(ANALYTICAL ANALYSIS) 



A. Low voltages 
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Substituting Ao given by Eq. (9) to Eq. (6), we arrive at 
the following equation for T: 



In particular, at V 



0, Eqs. (8) and (10) naturally 
yield Ao = A'o'^ and T = Tq. At sufficiently low voltages 
when {ri^'^''"'^eV/hujQ) ^ 1, the solutions of Eqs. (9) and 
(10) can be found analytically. In this case, Ao - 
Considering this, at low voltages, Eq. (10) yields 
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(13) 



where 77^'^ is the value of 770 at <C hujo/e , i.e., r/g'^ 
5 [IB]- As follows from Eq. (13), an increase in the applied 
voltage V leads to a decrease in the effective temperature 
of the electron- hole plasma (its cooling). 

Using Eqs. (5) and (13) at low voltages, we also obtain 
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(10) 



Here the current J is normalized by its characteristic 
value Jo = 2eLSo/T^"*"'' = 2eLif""'^ 



5 



B. Special cases 

In the special case V = fuvo/e (i.e., V ~ 0.2 V), from 
Eqs. (9), (10), and (12) we obtain 



T = 



In 



1 + % 



One can see that in this case 

1 



i.e., e = T and 



exp(^o/r) - 1 ' 



J = Jo 



(15) 

(16) 
(17) 

(18) 



At To = 300 K and 2L = 20 /^m, one obtains Jo 
32 A/cm. As follows from Eqs. (15) - (17) at V = hujo/e, 
T and Q tend to Tq if rj^'"^^ tends to zero, and T = 6, 
and they both increase proportionally to r\^'^°'^ as t/q^'^"^ 
tends to infinity. Thus, at V = Swo/e and r]Q^'^'^^ ^ 1, the 
effective temperatures are fairly high: T = 6 ~ hwo/ In 2 
(T = e ~ 3300 K). It is worth noting that J at F = 
fvjJo/e is independent of parameter r\^'^°''^ . 

In interesting (but nonrealistic) limiting case r\^'^°'^ = 
0, from Eq. (8) we immediately obtain A/q = AAq*. In 
such a case, T = Tq both at F = and V = tuvo/e. If 
V = hujQ/2e, one obtains 



To 



2[1 - (To/fiwo) ln(2 + 770) 



J - Jo(??o + l)exp( ) < Jo. 



(19) 



(20) 



At To = 300 K , Eqs. (19) and (20) yield, T ~ O.6OT0 = 
180 K and J ~ 31 mA/cm at Tq = 300 K and T ~ 
0.569To = 114 K and 0.64 mA/cm at Tq = 200 K. 



C. Long optical phonon decay time 

In the case of relatively long optical decay time when 
^decay ^ I y < ^g/e, neglecting terms of the order 
of A/;^Vj?^f^"''^ ~ exp(-;^a;o/To)/r?^"^"^ from Eq. (10) we 
obtain 



T~ 



eV 



In 



T^^o^'^^eV/hiuo) + miieV/HuJo) - 1] 



(21) 



At y = hujo/e, Eq. (21) yields the same value of T as 
Eq. (15) provided rj^^^^^ > M^". 
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FIG. 3. Electron-hole effective temperature T as function of 
applied voltage V for different values of r]^'^'^^ and To. Dashed 
lines correspond to T = 300K (upper panel) and T = 200 K 
(lower panel). 



V. EFFECTIVE TEMPERATURES AND 
CURRENT- VOLTAGE CHARACTERISTICS 
(NUMERICAL RESULTS) 



To obtain T—V and 8 — F dependences in wide ranges 
of parameter r\^'^°''^ and the applied voltage Eqs. (9) 
- (12) were solved numerically. Figures 3 and 4 show the 
voltage dependences of effective temperatures T and 
calculated for different values -q^^^'^'' and Tq = 300 and 
200 K. One can see from Fig. 3 that the electron-hole ef- 
fective temperature markedly decreases with increasing 
voltage, so that T < Tq in a certain voltage range [see 
also Eq. (13)], and then starts to rise. In the range of rel- 
atively high voltages {V < hwo/e), the T — V dependence 
is steeply rising [in line with analytical formula given by 
Eq. (21)] with T > Tq or even T » Tq. However it is 
intriguing that in a rather narrow voltage range where V 
is about some value Vd, the T—V dependences are am- 
biguous, so that these dependences as a whole are of the 
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0.1 

Voltage, V (V) 

FIG. 4. Optical phonon effective temperature Q as function 
of applied voltage V for different values of 7]^'^'^'^^ . Dashed line 
corresponds to = 300 K. 



S-shape. The appearance of the S-shape characteristics 
can be attributed to a decrease in parameter r]o with in- 
creasing T [sec Eq. (11)]. This corresponds to a decrease 
in To"*'"", and, hence, to an essential intensification of 
the intraband transitions, particularly, those associated 
with the rcabsorption of nonequilibrium optical phonons 
when T increases. This is because at high electron-hole 
effective temperatures the intraband transitions assisted 
by optical phonons take place between relatively high en- 
ergy states with their elevated density. When V exceeds 
some "disruption" voltage Vd, the net power acquired by 
the electron-hole plasma can be compensated by the in- 
traband energy relaxation on optical phonons only at suf- 
ficiently high T. As a result, in this case the clcctron-holc 
temperature jumps to the values corresponding to higher 
branch of the T—V dependence. Thus, the "observable" 
T — V dependences and their consequences can as usual 
exhibit hysteresis instead of the S-bchavior. One needs 
to point out that if the above temperature-dependent pa- 
rameter rjo is replaced in calculations by a constant, the 
calculated T — V dependences become unambiguous, al- 
though they exhibit a steep increase in the range V Vd. 

The behavior of T as a function of V markedly depends 
on parameters rj^'^'^"'^ and rjo. The width of the voltage 



range where T < Tq increases when parameter ^ be- 
comes smaller with increasing voltage. Simultaneously, 
the depth of the T—V dependence sag with T < Tq 
increases with decreasing tJq^'^"'^ as well as with decreas- 
ing Tq. At small t/q^ , the electron- hole cooling can 
be rather strong, particularly when Tq = 200 K. This is 
natural because faster decay of optical phonons prevents 
their accumulation (heating) and promotes the electron- 
hole plasma cooling when the Joule power is smaller than 
the power transferred from electrons and holes to optical 
phonons. It worth noting that the voltage range where 
the T—V dependence is ambiguous widens with increas- 
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FIG. 5. Normalized current- voltage characteristics for for dif- 
ferent values of r?*'^"'^- 
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As seen from Fig. 4, the optical phonon effective tem- 
perature also exhibits a S'-shape voltage dependence. 
However, contrary to the electron-hole effective temper- 
ature, > To at all the voltages under consideration. 
The values of Q at relatively high voltages steeply in- 
crease with increasing parameter Tjg*^""^. 

Comparing the T—V dependences calculated for differ- 
ent lattice temperatures, one can find that at moderate 
and large values of 770*^'"^, these dependences are virtu- 
ally independent of Tq. This is because in such a case 
the number of optical phonons A^ » Mq'^ and, hence, 
9 » To even at not too high voltages, so that the role of 
equilibrium optical phonons is weak. 

Invoking Eq. (5), the T — V dependences obtained 
above can be used to find the current- voltage characteris- 
tics. Figure 5 shows the J — V characteristics calculated 
using Eq. (5) and the T — V and Q — V dependences 
obtained numerically. As a consequence of the S-shape 
T — V and Q — V dependences, the J — V characteris- 
tics (as well as the voltage dependences of the dynamic 
conductivity considered in the following) are also of the 
S'-shape. According to Figs. 3 and 5, the T — V and 
J — V characteristics in the range of low and moderate 
voltages arc independent of parameter 7]^'^'^°'^ . This is 
in line with the results of the previous analytical analysis 
[see Eqs. (13) and (14)]. However, at relatively high volt- 
ages, distinctions in the J — V characteristics for different 
''-^'^"-y is significant although all of them tend to Jo when 



V approaches to Hujo/e. 



VI. DYNAMIC CONDUCTIVITY 

Knowing the T — V dependences, one can calculate 
the dynamic conductivity, cr^j, of a GL under the injec- 
tion pumping as a function of the signal frequency lu and 
the applied voltage V. To achieve lasing at the frequency 
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FIG. 6. Dependence of normalized dynamic conductivity 
Re (Ju)/(^o on signal frequency uj/2-k for different values of 
^decay — 3 X 10^^ s~^ (upper panel) and uq ~ 

15 X 10^^ s~^ (lower panel). 
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w, the real part of the complex dynamic conductivity at 
this frequency should be negative: Re cr^ < 0. As shown 
previously (see, for instance, Refs. [1,[3^), the interband 
contribution of the nonequilibrium electron-hole plasma 
with the quasi-Fermi energy fi and the effective temper- 
ature T is proportional to tanh[(?iaj — 2fi)/4T]. The in- 
traband contribution to Re , which corresponds to the 
Drude absorption, depends on uj, fi, and T as weU. It 
also depends on the time of electron and hole momen- 
tum relaxation on impurities and phonons r. The latter 
is a function of the energy of electrons and holes e. The 
main reason for the t — e dependence is a linear increase 
in the density of state in GLs with increasing e. In this 
case, = vaie/To), where i/q is the collision frequency 
of electrons and holes in equilibrium at T = Tq. Consid- 
ering this and taking into account Eq. (1), we can arrive 
at the following formula [l6| approximately valid in the 
frequency range a; 3> I'd : 



^^tanhf^^^Vc(-^^^ + ^f^^/^). (22) 

Here cto = /Ah and C = 2hi'o/TrTo. In high quality 
MGLs (with i/Q-i ~ 20 ps at To = 50 K [3,[3), assuming 
that vq oc To [s^, for To = 300 K one can set j/q — 
3 X 10^^ and hence, C ~ 0.005. For substantially 
less perfect GLs with fo ~ 15 x 10^^ s^^, one obtains 
C ~ 0.025. These data are used for the calculations of 
Re a^^. 

Figure 6 demonstrates the frequency dependences of 
Re (T^/ao calculated using Eq. (22) for different values 
of parameter rj^'^'^"'^ at different voltages V for fo ~ 
3 x 10^^ s~^ and j^o ~ 15 x 10"^^ s~^. As seen from Fig. 6 
at the injection conditions under consideration, the char- 
acteristic conductivity ctq is much smaller than the dc 
conductivity in the i-section (Ti^\i^=o = CToo- It is also seen 
that even at relatively large values of parameter rj^'^'^'^^ , 
the dynamic conductivity can be negative in the THz 
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range of frequencies provided the applied voltage is not 
so strong to cause the electron-hole plasma and optical 
phonon system overheating and the ambiguity of the volt- 
age characteristics. This is confirmed by Fig. 7. Figure 7 
shows Re CTtj/cTo as a function of the applied voltage. As 
demonstrated, the range of the signal frequencies where 
Re (Ttj < markedly shrinks and the quatity | Re cr^ \ de- 
creases when either ttiq'^'^"'^ or vq increase. In particular, 

at large values of ^ the achievement of the nega- 

tive dynamic conductivity and THz lasing can be compli- 
cated (at the temperatures Tq ~ 300 K when the optical 
phonon recombination mechanism dominates). This is 
because when t]^'^"'^ increases, the quantity becomes 
small. As a result, the optical phonon system is over- 
heated starting from relatively low voltages that leads to 
an "early" overheating of the electron-hole plasma [sec 
Figs. (3) and (4)]. If the value of Re is insufhcient to 
overcome the losses of the THz modes propagating along 
the GL structure, the structures with MGL can be used. 
In this case, the net dynamic conductivity of the MGL 
structure is given by Re (j,^ x if, where K is the number 

of GLs Wm- 



VII. LIMITATIONS OF THE MODEL AND 
DISCUSSION 

High density of the electron-hole plasma in the i- 
section under the injection conditions promotes the 
quasi-neutrality of this section. The recombination does 
not significantly affect the uniform distribution of the 
electron and hole densities in the i-section (assumed 
above) if the recombination current J is much smaller 
than the maximum injection current J^a which can be 
provided by the p- and n-sections. This imposes the 
condition J <C 2Jm. The latter inequality is equivalent 
to the condition that the recombination length is longer 
than the i-section length 2L. Considering a strong de- 
generacy of the electron and hole components in the p- 
and n-sections, respectively, the quantity can be es- 
timated as 



J j7 



[m? - (m. - ^yf 



(23) 



At To = 300 K, 2L = 20 ^m, [i, ^ 0.3 eV, and V = 0.1 V, 
one obtains Jq — 32 A/cm and 2Jm ~ 65 A/cm. As 
seen from Fig. 5, J <C Jq at least in the most interesting 
voltage range, where Re CTi^ < (see Fig. 7), the condition 
J ^ Jra is satisfied. 

In the above consideration we disregarded the leak- 
age current from the i-section to the p- and n-sections. 
This current includes the tunneling and thermionic com- 
ponents. Both these components depend on the height, 
^i, of the barriers between the p- and i-sections and i- 
and n-sections and the applied voltage V (see Fig. 2). As 
shown previously [37j , the tunneling current decreases 
with increasing width of the p-i- or i-n junction I [see 
Fig. 2(a)], because it is sensitive to the electric field at 



the junction [18|. The width in question depends on the 
geometrical parameters of the structure, in particular on 
the spatial distributions of donor and acceptors near the 
junction, and the thickness of the gate layer VFg, as well 
as the shape of the gates (in the structures with the elec- 
trical doping). So one can assume that this width can 
be sufficiently large to provide smooth potential distri- 
butions at the junctions. The effective height of the bar- 
riers at the p-i- and i-n-junctions, which determines the 
thermionic electron and hole current over these barriers, 
is equal to A = /ij — 2/i = [li — eV. It can be small 
at elevated values of the quasi- Fermi energy fj. which are 
necessary to achieve the negative dynamic conductivity 
(see, for instance, Refs. [1,13, |^, [l3)j i-e., in the most 
interesting case. This, in turn, implies that the electric 
field at the junctions and, hence, the tunneling current 
is decreased. In this case, the thermionic leakage current 
dominates over the tunneling leakage current, and the 
latter is disregarded in the following estimates. 

Taking into account the height of the barrier A, for the 
contribution of thermionic current through the junctions 
to the net terminal current one can obtain 



exp 



eV — fii 
T 



1 



T 

2evth^o{ 7^ 



exp( I - 1 



(24) 



Here Vth = {^2/tt^)vp. Comparing J and Jth given by 
Eqs. (5) and (24), respectively, one can conclude that the 
leakage (thermionic) current is small in comparison with 
the current associated with the optical phonon recombi- 
nation if [vthTo/ L) exp{—^i/T) <C exp(— fiwo/T) or 



hujQ > Tin 



L 



(25) 



Assuming that tq ^ 10^^^ s, factor (vthTo/L) in Eq. (21) 
is small when the length of the i-section 2L > 1 /.tm, i.e., 
at fairly practical values of 2L. Therefore, neglect of the 
thermionic leakage in our calculations in the above sec- 
tions is justified when /i^ > huQ ~ 0.2 eV (more precisely 
when Hi — htoo > T). At T > Tq, the latter inequality 
means that the donor and acceptor density in the perti- 
nent sections should be 



1 f UJQ 



3.3 X lO^^cm"^^ 



(26) 



In the case of the devices with the electrically-induced 
p- and n-section, the analogous condition sounds as 



Vg> 



4eWg f Wo 



ae \vf 



(27) 



Setting se = 4: and Wg = 10 nm, the latter condition 
corresponds to Vg > 1.5 V. 
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At a strong heating of the electron-hole plasma, say, 
a.t V = Hluq/c, using Eq. (15), condition (25) is replaced 

by 



\n{vthTo/L) 



In 



decay 

Vo 



eq 



decay 



(28) 



The latter inequality can impose somewhat stricter lim- 
itation on the values of /i, and than those given by 
Eqs. (25) and (26) if 77^"""^ > 1. 

It is notable that an effective confinement of the in- 
jected electrons and holes in the i-section by the barriers 
at p-i- and i-n-junctions in the GL structures under con- 
sideration can be realized by relatively low doping levels 
in p- and n-sections in comparison with structures with 
two-dimensional (2D) electron plasma in quantum wells 
on the base of the standard semiconductors. Indeed, the 
barrier height at = in the doped section of GL is 
equal to A = /i^ — eV. In standard 2D systems with the 
electron effective mass m* and the energy gap Ag, the 
barrier height is equal to A* = Ag -I- /i* — eV* , where 
/i* ~ nh^Tj* / m* . To achieve the same value of the quasi- 
Fermi energy /i in the standard 2D system as in GLs, one 
needs to apply the voltage V* =V + Ag/e. To provide, 
for example, the values — ^* = 0.3 eV, one needs 
Ei ~ 7.33 X 10^2 and S* ~ 3.82 x 10^^ ^^-2^ 

spectively. Thus, in the standard 2D electron systems the 
doping should be five times higher (for m = 4 x 10~^^ g). 
This is due to lower density of states in massless GLs 
near the Dirac point compared to that in the standard 
2D structures with m 7^ 0. Since the effective mass of 
holes M* in the standard semiconductors is markedly 
larger than to*, the realization of the barrier height at 
the p-i-junction sufficient for the effective confinement of 
electrons at elevated temperatures requires fairly heavy 
doping. Thus, in contrast to the standard 2D structures, 
the thermionic leakage current in the GL or MGL p-i-n 
structures under consideration can be sufficiently small 
without the employment of wide-gap p- and n-sections. 
In passing it should be mentioned that an extra confine- 
ment of the injected electrons and holes can be achieved 
if the p- and n-sections constitute arrays of graphene 
nanowires (doped or with electrically-induced high elec- 
tron an hole densities), so that the p-i-n structures con- 
sidered above are replaced by the P-i-N structures. 

One needs to stress that the assumptions (used in the 
above model) that the electron-hole plasma in the active 
region is virtually uniform as well as that the recom- 
bination current exceeds the leakage current are rather 
common in simplified models (the so-called rate-equation 
models) of in standard injection laser structures with the 
double injection (see, for instance, Refs. ^sl Isoj). 

As follows from the above calculations, the effective 
electron-hole and optical phonon temperature can be 
very high at V fkoo/e and be accompanied by ef- 
fects associated with the S'-shape characteristics. In this 



case, an expression for the rate of optical phonon decay 
(A/q — A/'Q'')/Tg^'^°^ used in Eq. (7) might be oversimpli- 
fied due to a strong anharmonism of the lattice vibration. 
Possibly, the effects related to a strong anharmonism can 
be taken into account by a proper choice (renormaliza- 

, • \ r , decay i decay t ,-\ i i i 

tion) ot parameters Tg and jjq . in the above treat- 
ment, for simplicity only one type of optical phonons with 
tiujQ ~ 0.2 eV was taken in to account. However, due to 
closeness of the optical phonon frequencies of different 
type in GLs and MGLs, the pertinent generalization of 
the model, adding computational complexity, should not 
lead to a marked change in the obtained results. Apart 
from this, at large effective temperatures, the radiative 
recombination and cooling (due to the radiative trans- 
fer of the energy outside the structure) can become es- 
sential [40l - |4^ . resulting in a limitation of these tem- 
peratures and affecting the S'-shape dependences. This 
means that considering the range or relatively high ap- 
plied voltages and, hence, strong injection, our purely 
"optical phonon" model should be generalized. How- 
ever, this concerns not particularly interesting situations 
in which the dynamic conductivity is not negative. 

As demonstrated, the main potential obstacles in the 
realization of negative dynamic conductivity and THz 
lasing in the injection GL and MGL structures at the 
room (or slightly lower) temperatures might be the intra- 
band photon (Drude) absorption and the optical phonon 
heating. These effects are characterized by parameters 
C cx i/Q and 770*^°°^ oc Tq'^'^"'^ , respectively. As for pa- 
rameter C, it can be sufficiently small in perfect MGL 
structures like those studied in Ref. , so the problem of 
intraband absorption can be overcome. However, if real 
values of parameter t]^^'^"'^ can not be decreased to an 
appropriate level (r\^'^°'^ ^1)) the achievement of room 
temperature THz lasing in the structures under consid- 
eration might meet problems. In the case of such a sce- 
nario, the utilization of lower temperatures, at which the 
recombination and energy relaxation is associated with 
different mechanisms, can become indispensable. 



VIII. CONCLUSIONS 

In conclusion, we have studied theoretically the effect 
of population inversion associated with the electron and 
hole injection in GL and MGL p-i-n structures at the 
room and slightly lower temperatures when the interac- 
tion with optical phonons is the main mechanism of the 
recombination and energy relaxation. In the framework 
of the developed model, the electron-hole and optical 
phonon effective temperatures and the current-voltage 
characteristics have been calculated as functions of the 
applied voltage and the structure parameters. It has been 
demonstrated that the injection can lead to cooling of the 
injected electron-hole plasma in the device i-section to 
the temperatures lower than the lattice temperature at 
low and moderate voltages, whereas the voltage depen- 
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dences can be ambiguous exhibiting the S-shape behavior 
at elevated vohages. The variations of the electron-hole 
effective temperature with increasing applied voltage are 
accompanied with an increase in the optical phonon ef- 
fective temperature. Using the obtained voltage depen- 
dences, we have calculated the dynamic conductivity and 
estimated the ranges parameters and signal THz frequen- 
cies where this conductivity is negative. The electron- 
hole cooling might substantially promote the realization 
of THz lasing at elevated ambient temperatures. In sum- 
mary, we believe that the obtained results instill con- 
fidence in the future of graphene-based injection THz 



lasers although their realization might require a thorough 
optimization. 
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